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Effects of cyclic fluctuations in meridional circulation using a low 
order dynamo model 
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ABSTRACT 

We develop and subsequently explore the solution space of a simple flux transport dy- 
namo model that incorporates a time dependent large scale meridional circulation. Based on 
recent observations we prescribed an analytical form for the amplitude of this circulation and 
study its impact in the evolution of the magnetic field. We find that cyclic variations in the 
amplitude and frequency of the meridional flow affect the strength of the solar cycle. Varia- 
tions in the amplitude of the fluctuations influence the shape of the solar cycle but are only 
relevant to the cycle's strength variations when they occur at a frequency different from or out 
of phase of the solar cycle's. 
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1 INTRODUCTION 

Our civilization is increasingly becoming more and more depen- 
dent on energy distribution, communication networks and satellite 
operation since these technologies provide an important backbone 
for our daily activities. Nevertheless, these key technological as- 
sets may face dangers that one would like to minimize. Our Sun, 
so important for the life on our planet, is now being pointed as a 
source of problems for these technologies. Solar magnetic storms 
can surreptitiously hit Earth and damage all the structures previ- 
ously mentioned. At the origin of these storms we can find the 
large scale solar magnetic field. This field is believed to be origi- 
nated by a magnetohydrodynamic dynamo that converts kin etic en- 
ergy from the solar plasma motions into magnetic energy jParker I 
Il955l) . Several dynamo models have been developed t o explain the 
main observational feature s of the solar magnetic field djouve etal\ 
120081 : ICharbonnead feoiO). These models can reproduce the mag- 
netic field polarity reversals (every 1 1 years) and many field spa- 
tial features. Touted as being the most promising of the several 
existing types, flux transport dynamo models have been tenta- 
tively use d for the first time as a tool for prediction s of the so- 
lar cycle dDikpatief fl/]|2006t IChoudhuri erq/ll2007l) . These pre- 
dictions were the firsts to use full dynamo models to forecast the 
behavior of the solar cycle. The models solve the mean-field ax- 
isymmetric equations for the magnetic field evolution on a back- 
ground structure that incorporates parameterized physical mech- 
anisms such as solar rotation, magnetic diffusivity, etc. The de- 
nomination "flux transport models" comes from the fact that they 
incorporate the solar meridional circulation, a conveyor belt-like 
plasma flow that carries magnetic field from the equator to the 



poles near the surface and from the po les towards the equator in 
the base of the solar convection zone dDikpati and Charbonneau I 
ll999l;IChatteriee et aZl200llMunoz -Jaramillo et a 320091) . The am- 
plitude or strength of this circulation c ontrols the period and ampli- 
tude of the produc ed magnetic field dNandv and Choudhuri 112002 



Lopes and Passos 20 09a|). With ju st a few exceptions dRempel I 



20061 . Karak and Choudhuri 1 1201 lh . most of these models work 
on the kinematic regime, which implies that the amplitude of 
the meridional circulation is not affected by the electromagnetic 
Lorentz force feedback of the magnetic field on the flow. Obser- 
vational evidence of meridional circulation are very hard to ob- 
tain and current values for the surface poleward average veloc- 
ity is centered around 10 to 20 ms _1 . The most reliable data is 
available onl y for the last couple of decades an d the recent mea- 
surements of lHathawav and Rightmirel d2010al) indicate that the 
strength of this flow might have changed by about 25% from the be- 
ginning of cycle 23 to cycle 24. Experimental evidenc es for a vari- 
able m eridional circulation were already reported by (Ko mm et al\ 
I1993I) . Also, a recent inversion methodology based in a simpli- 
fie d dynamo model and the annual sunspot time series proposed 
bv lPassos and Lopes ld200ll ): |Passos (2012a) suggests that the am- 
plitude of the meridional circulation changes significantly from cy- 
cle to cycle and that those variations could explain (partially) the 
observed solar variability. This information is not usually taken 
into account in dynamo based predictions but its relevance is 
now rec ognized and is starting to be addressed by some research 
groups (Karak1 l20ld iHotta and Yokoyama 2010; Nandye/a/ 
1201 lb It ha s been shown bvlMininni et al\ d2000h : lwilmot-Smith 
(2005) and IPassos and Lopes I ((2008, 2011) that a truncated ver- 
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sion of the flux transport dynamo equations set or low order dy- 
namo models can be used as a first order approximation to study 
the temporal behavior of the solar magnetic field. These ID trun- 
cated models allow to calculate the evolution of the magnetic field 
strength by taking into account the main physical mechanisms in a 
simplified way in the kinematic regime, without the need to incur in 
heavy numerical calculations. Exploratory studies made with these 
models are very useful for identifying and to help focusing in the 
relevant aspects and physical mechanisms that should be studied in 
depth by 2D numerical mo dels. In this work we bu ild upon the low 
order dynamo presented in lPassos and Lopes by taking into 

account a time dependent meridional circulation's amplitude pro- 
file. After the derivation of the new low order model equations, we 
p arameterize t he me ridional circulation based on the observations 
of lHathawav I d20 1 lh and study the impact of this time dependence 
in this dynamical system's solution. The final section is dedicated 
to comments and remarks about the results. 



2 THE MODEL 

As previously mentioned, the model presented here is a vari- 
at ion of the kinematic lo w order dynamo model de veloped 
in lPassos and Lopesl j2008h and[Passos an d Lopesl d201ll) . In this 
version we consider that the amplitude of the meridional circulation 
is time dependent and the derivation of the low order model fol- 
lows this directive. We start with the equatio ns for a flux transpor t 
mean field axisymmetric dynamo as shown in lCharbonneaul J2010h . 
These equations give us the evolution of the mean solar magnetic 
field, B = + B,,, classically decomposed into its toroidal, B^, 
and poloidal, B /; = V x (A p e^,) components with A p representing a 
potential vector field. 
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where we have r = r sin 9, VC1 represents the differential rotation of 
the Sun, v ; , is the flow in the meridional plane and r] is the magnetic 
turbulent diffusivity. One of the simplifications used in this model is 
to assume an average magnetic diffusivity for the entire convection 
zone (drj/dr = 0) and plasma incompressi bilit y. Following the sug- 
gestions found in lKitchatinov et a/I J2000h and lPontieri et a\ b003h 
we add an extra term, T ~ yB^/Snp, to account for magnetic flux 
removal by magnetic buoyancy. Here y is a constant related to the 
buoyancy regime and p is the plasma density. As usual in these 
models the regeneration mechanism from toroidal to poloidal field, 
the so called tt-effect is represented by a. In this case, for sim- 
plicity, we do not consider any non-linearity in a. In the following 
steps we assume that \ p depends ex plicitly on time and w e use the 
dimensional approach suggested by iMininni et al 1 boooh to trun- 
cate the dynamo equations by substituting V — > I/Cq, where t a is a 
specific length of interaction for the magnetic fields, usually taken 
as £ ~ 0AR o . This truncation ensures that we are bounding our 
solution space to magnetic phenomena that occur in the scale of £ , 
presumably the large scale solar magnetic field responsible for the 



solar cycle. After grouping terms in and A p we get 
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These coefficients now contain all the structure parameters in the 
model, i.e. c\, c 2 and c 3 assume the role of magnetic diffusivity, ro- 
tation and buoyancy respectively in this low order dynamo model. 
Next, we take the derivative of equation ^ with respect to time and 
drop the A p dependence by substituting equation l[4j in the terms 

dA 

with and by noting that cjAp can be extracted from equation 
<[3j . After this mathematical workout we finally get that 
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with a + 0. The solution's space of this dynamical system is de- 
fined by the structural coefficients c„ and by the analytical form of 
v p (t). As a note, it is important to refer that in this kind of reduced 
systems, the units in which some of the quantities are presented do 
not always coincide with the real units, rendering the direct appli- 
cation of known physical values troublesome. Nevertheless studies 
based on the relative variation of these parameters can be done and 
this work follows that line of thought. 

A "static" \>p reference solution for equation {8} is calculated 
by assuming a constant meridional flow amplitude. The values used 
for the coefficients c„ in this static solution were found by fitting 
e quation (QQ to the construct ed proxy for the toroidal field presented 
in Pass os and Lopes (for further details about the methodol- 



ogy used you can also see lLopes and Passos (2009a)). Using these 
c„ values the system presents a solar-like solution with S^, a proxy 
representation of the solar cycle, showing a cyclic behavior with a 
period of about 1 1 years (figure (TJ). In the parameter regime used 
in this work, v p (t) behaves mathematically as a source term and ci, 
the diffusivity as a sink term. The former is one order of magnitude 
higher than the latter. Also, in this parameter regime, magnetic flux 
removal by buoyancy, c 3 , is the main saturation mechanism in place 
to avoid field growth. In the following of this work we bound our- 
selves to the study of variations in v p (t) maintaining the other coef- 
ficients with the values presented in the static solution. A complete 
study of the full parameters space is therefor deferred to a future 
work but some preliminary results are presented in the conclusions 
section. 



2.1 Introducing cyclic fluctuations in v p (t) 

According to the latest measurements of the surface merid- 
ional flow amplitude spanning a full magnetic cycle from 
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Figure 1. Static solution (obtained for constant v p ), after the system evolved 
for some time into stability. The solid black line represents B?, our chosen 
proxy to represent the solar cycle and the dashed blue line is the scaled 
amplitude of v p /Cq (blue) in order to plot both quantities in the same scale 
(in this example the scaling used is 1000 |v„|). Values used here C[ = -0.01, 
c 2 a = -0.095, c 3 = 0.002, v p /£ =-0.l. 

lHathawav and Rightmire (2010a), v p (t) varies in a roughly sinu- 
soidal way, reaching a maximum amplitude near the half of the de- 
creasing part of the solar cycle and droppin g to its lower value near 
the sun spot maximum (see figure (4) of lHathawav and Rightmire] 
Based in this result we propose an Ansatz where we take 

oc v p o+ A sm(wt+8) where A is the amplitude of the meridional 
flow fluctuations, a> is the frequency of these fluctuations and 9 is 
used to control the initial phase. This parameterization translates in 
a fluctuation of v p around a mean value of v p0 (value used in the 
static solution). The top panel of figure $2% shows the solution for 
oj = ~ with T = 11 years and a fluctuation's amplitude A=0.025 
which corresponds to an amplitude variation of 25% v p0 (roughly 
the amplitude variation bet ween cycle minimum and maximum for 
solar cycle 23 presented in Hat haway and R ightmire feOlOal) ). It 
is important to mention that even if our meridional circulation is 
time dependent, the model is still operating in the kinematic regime 
since v p does not depend on the magnetic field. 

We notice that the initial phase, 6, has no impact on the so- 
lution's shape (for the range of tested parameters). One interest- 
ing result is the fact that for a> = M the phase difference between 
v p and B^ is the approximately the same as the observed one, i.e., 
maximum amplitude of v p at the decreasing phase of the solar cy- 
cle and minimum amplitude of v„ near the cycle maximum. We 
used as initial conditions for solving equation that S^(0)=0.01 
and dB,f,/dt(Q)=0. For these specific initial conditions, B^t) enters 
a steady oscillation regime approximately after t= 100. For an os- 
cillation amplitude of v p of 25% and an initial phase = 0, the 
phase difference between B^ii) and v p "locks" around t=300 while 
for 6 = n/2 this occurs at approx. t=750. For different initial condi- 
tions we get different times for the "phase lock". For higher values 
of v p (either v p o or A) the time to achieve the phase lock decreases. 
The phase lock occurs when the frequency a> associated with the v p 
fluctuation is the same as (or very close to) the natural frequency 
of B$. In this low order model, the period of the cycle is given pri- 
marily by c-ict w ith a small influence of c\ - v p /£o (for details c.f. 
IPassos I J2012al) ). This small dependence on v p seems enough to 
ensure that after some time fi^ is synchronized with v p {t). Another 
way of explaining this is to resort to a {v p , B$] phase space. If the 
solution in this phase space is a limit cycle then the two quanti- 
ties will synchronize phases after the solution evolves towards the 
attractor. See figure ([3} for an illustrative example. 

With this set of parameters the system is well behaved and 
has a stable solution in the form of an attractor or limit cycle, best 
viewed in the phase space of B(t) (figure (O bottom panel). On the 
other hand different values of the fluctuation's frequency, a>, yields 
some impact in the solution. For fluctuations with frequencies dif- 
ferent than that of the solar cycle, like the one presented in the 
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Figure 2. The lines representation is the same as in figure (1). Top: c\ = 
-0.01, c 2 a = -0.095, c 3 = 0.002, v p0 =-0.1, A=0.0225, (O = 2n/U and 
= 0. Middle: Same parameters as top panel, but with a change of the 
fluctuation's period, a) = 2n/\4. In the bottom panel are represented the 
phase space B<p, dB^/dt for both solution with T=l 1 (left) and T=14 (right) 
sampled at regular intervals. 




Figure 3. [v p , B<»| phase space for the solution between t = 1000 and 
t = 1 100 using w = 1 1 (left) and a> = 12 (right). The limit cycle on the left 
panel is observed because at and the frequency of oscillation of Bj, is the 
same. 

middle panel of figure (O, we observe a clear modulation of the 
field strength. This is best viewed in the corresponding phase space 
where we change from a stable well defined limit cycle to a "limit 
region". A natural variability appears in the system even if the av- 
erage meridional circulation amplitude remains constant (here v^o)- 
Solutions computed with different values of the fluctuation's am- 
plitude, show that it has an influence mainly in the shape of the 
cycle, creating higher asymmetries between its rising and falling 
parts. We also observe very small changes in the frequency of the 
cycle but there are no signs of any long term variability (ampli- 
tude variations in the cycle's strength). With fluctuation's as high 
as 200% v p o, the solution in the phase space remains a limit cycle. 
According to observations, the most realistic scenario is to con- 
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Figure 4. Solution assuming a variation in the fluctuation's amplitude A, 
and frequency, to at a certain moment. Fixed parameters c\ = —0.01, ciot = 
-0.095, c 3 = 0.002, 6 = 0, v,, =-0.1. In the top panel we used oj = 2tt/11 
and A=0.025 (black) until t = 1052 and cj = 2tt/13 and A=0.05 (gray) 
afterwards. The correspondent phase space is presented in the bottom panel. 
The line in black corresponds to the solution between t = 1000 and t = 1052 
and the gray dots are a sampling of the solution at regular intervals after the 
change. 

sider small variations in the fluctuation's amplitude and perhaps in 
the period (of the order of a couple of years). These variations will 
create a variability in the strength difference between cycles N and 
following cycles. In this case the observed solar cycle variability 
becomes dependent of changes in the fluctuations pattern of v p . We 
present in figure © a test case where the amplitude of the fluctu- 
ations rise from 25% v p o to 50% v,,o, and the frequency decreases 
from u) = 2^/11 to to = 2n/l3 at a t=1055. In this case we find 
that both amplitude and frequency of the fluctuations influence the 
system's response. Changing the oscillation amplitude only (main- 
taining the frequency) results in a small variation of the limit cy- 
cle, i.e., a very smooth variability in the amplitude of the cycle 
appears (~0.2%) and phase lock between flow and field is main- 
tained. If one allows for variations in the oscillation frequency as 
well, a larger variability of the order of 30% or larger appears (tran- 
sition from limit cycle to limit region in the phase space). When the 
fluctuation's frequency occur out of phase with the solar cycle, the 
amplitude of the fluctuations influences the strength of the future 
cycles in a more pronounced way. 



3 CONCLUSIONS AND FINAL REMARKS 

The main objective of this work is to quickly probe the impact that 
cyclic fluctuations in the solar meridional circulation profile have 
in the dynamo process operating in our Sun. To to so, we use a 
simplified dynamo model where meridional circulation amplitude 
is forced in a sinusoidal way that mimic recent observations. We 
found that regular cyclic fluctuations in the amplitude of the solar 
meridional circulation do not seem to have an impact in the overall 
inner works of the solar dynamo. This is specially true when the 
period of these fluctuations is the same as the solar cycle's. Nev- 
ertheless, if the frequency at which the meridional flow varies be- 
comes different from the natural frequency of the solar cycle, then 
the cycle becomes naturally variable. In terms of phase space of the 



magnetic field we go from a well behaved solution, a limit cycle at- 
tractor, to a strange attractor (or attracting region). This might be 
a characteristic of this specific model because here the meridional 
flow term, v p acts as a source term. When this source term is not 
synchronized with the natural frequency of the solar cycle set by 
c 2 a, the variability appears. In this reduced model, the relationship 
between cycle period and amplitude that results from changes in 
the meridional flow presents a non typical solar behavior. Observa- 
tions show that, in average, the stronger solar cycles have shorter 
periods (amplitude-period rule) while in this model we get stronger 
cycles having longer periods. This apparent shortcoming can give 
us some clues about other quantities that can be varying over time. 
More specifically, variations in the physical mechanisms that are 
present in c 2 a, could produce the desired effect. Variations in the a 
effect or in the solar rotation could occur in parallel with meridional 
flow changes to produce the cycle's amplitude-period rule. 

According to this model, amplitude variations in the fluctua- 
tion profile of v p have an impact in the shape of the solar cycle, in- 
creasing the asymmetries between rising and falling parts and even 
creating double peaked cycles. In this case a solar-like feature is ob- 
served, i.e., stronger cycles tend to have a steeper rising phase than 
weaker cycles. Large and long term variations in the amplitude of 
the solar cycle only occur when the frequency and amplitude of the 
fluctuations change in parallel. 

Interesting to note is the fact that when our dynamo equation 
is forced with an sinusoidal meridional flow, the phase difference 
between both flow and field is nearly the same as the observational 
one. This is true as long as c 2 a is set to reproduce the 1 1 year peri- 
odicity of the solar cycle. We find that for different cia values (e.g. 
variations in the solar rotation or in the case of other stars rotating 
faster or slower than the Sun) the phase difference between merid- 
ional flow and magnetic field changes. As long as both frequencies 
(a> and the solar cycle's) are the same (or very close) the phase dif- 
ference between flow and field seems to lock independently of the 
initial phase used. Another parameter that influences the phase 
difference between field and flow is C\. By changing the diffusivity 
of the system the phase lock between field and flow can be modi- 
fied. 

Worth mention is the fact that the exact moment, in respect to 
cycle N, at which the variation in the fluctuations regime occurs has 
an impact in the amplitude of the following cycles. Depending on 
the chosen variation scenario (amplitude, period or both of v p ) cycle 
N+ 1 can be stronger or weaker than cycle N. A similar effect as also 
been r ecently reported in 2D dynamo simulations by Nan dv etal\ 
J201 J). We defer the details of this and other effects for a future 
work where we plan to perform a complete parameter space study 
of the model (including variations in c\ and C3). 

From a physical point of view the information that can be ex- 
tracted from such a truncated model is limited. Nevertheless we get 
important clues about the system's overall behavior when forced 
under certain parameterizations. One of the questions that this 
model in this present kinematic form does not address is what could 
be the cause(s) of the observed variations in the meridional flow. A 
possible explanation could be that the meridional circulation, be- 
ing a weak flow, can be influenced by the Lorentz feedback from 
the magnetic field. This feedback can be enough to modify sig- 
nificatively the meridional circulation. This scenario is supported 
by a recent analysis of the large-eddy glob al MHD simulat ions of 
the solar convection zone produced by Ghizaru et al 1l l2010l) where 
IPassos et al\ d2012b!) find evidence that the toroidal field at the base 
of the convection zone modulates the amplitude of the meridional 
circulation. If meridional circulation fluctuations are produced in 
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this way, they should occur at the same frequency as the solar cy- 
cle. Although different variations in amplitude from cycle to cycle 
could occur, the phase difference between this flow and the mag- 
netic field should remain the same. Future observation will gives us 
the answers. 

On the other hand, if indeed the Lorentz-feedback of the field 
into the flow is the reason (or partially responsible) for the ob- 
served amplitude fluctuations then it is reasonable to assume that 
this same feedback will also influence the solar rotation, although 
on a smaller scale. According to the found results, variations in 
these two large scale flows would be more then enough to produce a 
variable solar cycle. More details about this physical scenario could 
be presented here but those would be more speculative. Our inten- 
tion here is not to create/feed speculations but to explore plausible 
physical scenarios and motivate future studies. 

The final remark that we would like present is the fact that, 
since numerical dynamo models are now starting to emerge as fore- 
casting tools for solar activity, meridional variation mechanisms 
should be studied/implemented in order to improve their reliabil- 
ity. Most probably this will require a departure from the classical 
kinematic approach. In terms of the future, we also believe that if 
we keep monitoring v p (t) (and probably £l(t)) in the Sun then, ac- 
cording to the presented model, we would be able to predict the 
behavior of future solar cycles since the variability associated with 
the large scale flows seems to be mostly deterministic. 
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